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Abstract. The general normal ordering problem for boson strings is a combinatorial 
problem. In this talk we restrict ourselves to single-mode boson monomials. This 
problem leads to elegant generalisations of well-known combinatorial numbers, such 
as Bell and Stirling numbers. We explicitly give the generating functions for some 
classes of these numbers. Finally we show that a graphical representation of these 
combinatorial numbers leads to sets of model field theories, for which the graphs 
may be interpreted as Feynman diagrams corresponding to the bosons of the theory. 
The generating functions are the generators of the classes of Feynman diagrams. 
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1. Boson Normal Ordering 



In this note we give a brief review of the combinatorial properties 
associated with the normal ordering of bosons, and the model Feynman 
graphs which result. 

Combinatorial sequences appear naturally in the solution of the 
boson normal ordering problem PP, 

The normal ordering problem for canonical bosons [a, a'] = 1 is re- 
lated to certain combinatorial numbers S(n, k) called Stirling numbers 
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of the second kind through [3] 

n 

( a t a ) n = ^S(n,k)(at) k a k , (1) 
k=l 

with corresponding numbers B(n) = Y^l=i S( n > called Bell numbers. 
In fact, for physicists, these equations may be taken as the definitions of 
the Stirling and Bell numbers. For quons (q-bosons) satisfying [a, a'] q = 
aa^ — qa^a = 1 a natural q-generalisation _4j of these numbers is 

n 

(^a) n = Y,S q (n,k)(a^ k a k . (2) 

k=l 

In the canonical boson case, for integers n,r,s > we define generalized 
Stirling numbers of the second kind S rtS (n, k) through (r > s): 

ns 

[(a^YaT = (a^^ ]T S r , s (n, k){a)) k a\ (3) 

k=s 

as well as generalized Bell numbers B rtS (n) 

ns 

B r , s (n) = Y J SrA^ k )- ( 4 ) 

k=s 

For both S r ^ s (n,k) and B r ^ s (n) exact and explicit formulas have been 
found ^ E] • We refer the interested reader to these sources for further 
information on those extensions. However, in this note we shall only 
deal with the classical Bell and Stirling numbers, corresponding to 
B\ } i(n) and S\^{n) in our notation. 



2. Generating Functions 

In general, for combinatorial numbers g(n) we may define an exponen- 
tial generating function G{x) through 

OO Yl 

G{x) = Y J 9{n) X -. (5) 

n=0 

For the Bell numbers, the generating function takes the particularly 
nice form|Hj 

00 x n 

G(x) = Y Bin)— = exp(exp(x) - 1). (6) 

n=0 
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at a (a+af=a+a + a^a 2 (a + a) 3 =a + a + 3a +2 ef + a* 3 a 3 



B(1)=1 B(2)=2 B(3)=5 



Figure 1. Arrow graphs for (a^a) n n — 1,2, 3. 

Some initial terms of the sequence {B(n)} are {1, 2, 5, 15, 52, 203, 877, . . .}. 

3. Graphs 

A convenient way of representing combinatorial numbers is by means 
of graphs. To illustrate this, we now consider a graphical method for il- 
lustrating the combinatorial numbers associated with the normal order 
expansion of (a)a) n . 

A single arrow represents a "time-segment" of a line corresponding 
to the "propagator" (a^a). Thus we may concatenate one, two or more 
arrows to form a single line, or propagator (at a). However, two lines 
correspond to two distinct propagators a) 2 a 2 , and so on. Further, the 
constitutent arrows are labelled, for example by time, and so they may 
only be concatenated respecting the time ordering. These rules are 
illustrated by the diagrams of Figure 1, in which we consider the cases 
of 1, 2, and 3 arrows respectively. We have pre-emptively labelled these 
numbers as Bell numbers - which fact we demonstrate below. For the 
case of 4 arrows (Figure 2) we have additionally given the individual 
associated symmetry factors (in fact Stirling numbers) which add to 
B(4) = 15. It should be clear from this illustration how the time- 
ordering rules are applied to give the symmetry coefficients. Thus there 
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(a+af = st a + 4a* 2 a 2 + 3a +2 a* + Ba^a 3 



+ 



S(4,1)=1 



S(4,2)=7 



S(4,3)=6 



S(4,4)=1 



Figure 2. Arrow graphs for (a* a) 4 . 

is only one way in which we can concatenate 4 arrows respecting time- 
ordering (first grouping) , 4 ways in which we can divide the arrows into 
a set of 3 and 1, and so on. 

From these first few examples it would seem that these graphs are 
essentially like the Feynman Diagrams of a zero-dimensional (no inte- 
gration) Model Field Theory associated with H = a) a. In other words, 
at order n the total number of graphs would appear to be B(n) while 

the individual coefficients of a k are given by S(n,k). In order to 
show that this is indeed the case, we must be able to count the number 
of graphs associated with a given number n of arrows. To do this we 
can use the First of Three Great Results. 



This First Great Result is sometimes known as the connected graph 
theorem [Jj. It states that if C(x) = Y^n=i c(n)x n /n\ is the generating 
function of labelled connected graphs, viz. c(n) counts the number of 
connected graphs of order n, then 



4. First Great Result 




(7) 



is the generating function for all graphs. 
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We may apply this very simply to the case of the arrow graphs 
above. For each order n, the connected graphs consist of the single 
graph obtained by concatenating all the arrows into one propagator. 
Therefore for each n we have c(n) = 1; whence, C(x) = exp(x) — 1. 
It follows that the generating function for all the arrow graphs A(x) is 
given by 

A(x) = exp(exp(x) - 1) (8) 

which is the generating function for the Bell numbers. 

Such graphs may be generalised to give graphical representations for 
the extensions B r ^ s {n) jSJ. 

However, just as an abstract group is capable of more than one 
presentation, there are many graphical representations for a given com- 
binatorial sequence; and we now give an alternative one for the numbers 
B(n),S(n,k) due to Bender and collaborators jOj, [TUj . 



5. Second Great Result 

As before, we shall be counting lines. A line starts from a white dot, the 
origin, and ends at a black dot, the vertex. What we refer to as origin 
and vertex is, of course, arbitrary. At this point there are no other rules, 
although we are at liberty to impose further restrictions; a white dot 
may be the origin of 1,2,3,... lines, and a black dot the vertex for 1,2,3,... 
lines. We may further associate strengths V s with each vertex receiving 
s lines, and multipliers L m with a white dot which is the origin of m 
lines. Again {V s } and {L m } play symmetric roles; in this note we shall 
only consider cases where the L m are either or 1. 

We illustrate these rules for four different graphs corresponding to 
n = 4. 

There is an generating function G(x, V, L) which counts the number 
g{n) of graphs with n lines arising from the above rules |llj : 

oo m jm oo s 

G(x,V,L) = ex P (£L m -— ^£^1 

m=l " s=l y=0 

- E*(*)^ ( 9 ) 

n=0 

Consider the following example: L m = for all m ^ 1, that is, we 
allow only one line from each origin (with multiplier 1); there is no 
restriction on the number of lines to a vertex, and V s = 1 for all s. We 
give an example of the cases n = 1,2,3,4 in Figure 4. Note that for 
correct counting the lines should be labelled, as they were in the case 
of the arrows above. 
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L ,m o v > 




-o v , 
-o v , 



(L/V t V 3 (LfLJVf 



L 4 V 4 



0-/ (V/ 



Figure 3. Some examples of 4- line graphs. 



The generating function G(x) which counts the lines corresponding 
to the above rules follows immediately from Eq.@ 



00 x m d m 00 y s 

G(x) = exp(£L m — )ex P Q> s -) 

m=l a S=l y=Q 

1 oo s 

= exp(xd/dy)exp(e y - l)\ y=Q 

oo n 

= exp(e*-l) = J2B(n)^j. (10) 

n=0 



The penultimate step is a consequence of the Taylor expansion. 

We thus have yet another representation of the integer sequence 
{B(n)}. Note that when L = {L m } and V = {V s } are integer sequences 
we obtain an integer sequence from Eq.Q. A convenient method of 
obtaining the resulting integer sequence is afforded by the next useful 
result. 
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n= 1 



n = 2 



S(1,1) = 1 B(1)=1 



• o 

• o 



• o 

n = 3 m o 



S(2, 1) = 1 S(2,2) = 1 B(2) = 2 



-o 




S(3, 1) = 1 3(3,2) = 3 S(3,3) = 1 B(3) = 5 

Figure 4- Graphs of second type for B(n), n — 1,2, 3. 

6. Third Great Result 

Straightforward manipulation of series shows the following: Define 

CO Yl 

A i {x) = Y J ^i{n)— i = 1,2,3. (11) 



ra=0 



??.! 



Then if 
we have 



A 1 (xd/dy)A 2 (y)\ y=0 = A 3 (x) (12) 



a 3 (n) = ai(n)a 2 (n). (13) 



This is a useful and rather surprising equality. 

Using the results of the previous section, Ea. lfTBT) enables us to create 
graphical representations of products of integral sequences. For exam- 
ple: if in the case above we chose L m = 1 for all m, enabling any number 
of lines from each origin (with multiplier 1), the resulting sequence of 
graphs would have given us a representation of the integer sequence 

{B(n) 2 }M- 

We exemplify the possibilities offered by application of the Third 
Great Result by our two final examples. 

Example 1: Generating function for the sequence {B{n)B(n + 1)}. 
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From Eas. (|ll|) and the required generating function is given 

by 

G(x,V,L) = A l (xd/dy)A 2 (y)\ y=Q (14) 

where A\{x) is the generating function for {B(n)} and A 2 (x) is that 
of {B{n + 1)}. Note that A\(x) = exp(exp(x) — 1) from Eq.(jHJ), while 
A 2 (x) = {d/ dx)A\{x) = exp(exp(x) — 1 + x). 

This shows that the graphs for the sequence {B(n)B(n + 1)} may 
be obtained by putting L m = 1 for all m, so that there are any number 
of lines emanating from an origin, and with multiplicity 1. For the 
sequence {B(n + 1)} we have V s = 1 for all s ^ 2 and V 2 = 2, so that 
any number of lines may end at a vertex, and all have strength 1 except 
for the case where two lines meet at a vertex, when the strength is 2. 
The generating function is, according to Eq. ([13j) . 

OO n 

G 1 (x,V,L) = Y / B(n)B(n + l)^r. (15) 

It may be explicitly obtained after some formal algebraic manipulation 
based on the Dobihski formula ^5] 

1 00 h n 

as 

OO -I 

G 1 (x J V, J L) = ^T ? exp(exp((fc + l)x)-2). (17) 

k=0 

The formal series Q17|l diverges for all x > 0, although having finite 
Taylor coefficients for x = 0. Such formal series are nevertheless useful 
in representing integer sequences. 

Example 2: In our last example we retain all the derivative terms in 
Eq.(|ll|) but choose V\ = V 2 = 1; V s = 0, s > 2, thus allowing vertices 
where at most two lines meet. The corresponding generating function 
is defined by 

oo n 

G 2 (x,V,L) = Y,B(n)I(n)- (18) 

where the Involution numbers I(n) = 1,2,4,10,26,76,... are defined 
through their generating function 

2 oo n 

G 7 (x)=exp(x + ^) = ^/(n)^ (19) 
and are special values of the Hermite polynomials H n {x) 



I(n) = #„(-L)/(-V2i 
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Consequently 




= E^exp(Ml + ^)-l). (21) 

In obtaining Eq. (|21j) we have used the standard form of the generating 
function of the Hermite polynomials. Again, we must consider the series 
Eqs. ([20|) and (|21|) as formal power series, since for example they diverge 
for x > 0. 

In conclusion, we emphasize that the expressions of Eqs. (|17|) and 
(f2l|) constitute exact solutions of Model Field Theories defined by 
the appropriate sets {V s } and {L m }. Many other applications and 
extensions of the ideas sketched in this note will be found in [S]. 

Acknowledgements 

We thank Carl Bender and Itzhak Bars for interesting discussions. 



References 

1. Blasiak, P., Penson, K.A. and Solomon, A.I.: The general boson normal 
ordering problem, Phys. Lett. A 309 (2003), 198. 

2. Blasiak, P., Penson, K.A. and Solomon, A.I.: The boson normal ordering 
problem and generalized Bell numbers, Ann. Comb. 7 (2003), 127. 

3. Katriel, J.: Combinatorial aspects of boson algebra, Lett. Nuovo Cimento 10 
(1974), 565. 

4. Katriel, J.: Bell numbers and coherent states, Phys. Lett. A. 273 (2000), 159. 

5. Comtet, L.: Advanced Combinatorics, Reidel, Dordrecht, 1974. 

6. Wilf, H.S.: Generatingfunctionology, Academic Press, New York, 1994. 

7. Stanley, R.P: Enumerative Combinatorics vol.2, Cambridge University Press, 
1999; Riddel, R.J. and Uhlenbeck, G.E.: On the theory of the virial development 
of the equation of state of monoatomic gases, J. Chem. Phys. 21 (1953) 2056; 
Aldrovandi, R.: Special Matrices of Mathematical Physics, World Scientific, 
Singapore, 2001. 

8. Blasiak, P., Duchamp, G., Horzela, A., Penson, K.A. and Solomon, A.I.: Model 
combinatorial field theories, to be published. 

9. Bender, CM, Brody, D.C. and Meister, B.K.: Quantum field theory of parti- 
tions, J. Math. Phys. 40 (1999) 3239; Bender, CM., Brody, D.C, and Meister, 
B.K.: Combinatorics and field theory, Twistor Newsletter 45 (2000) 36. 

10. Bender, CM. and Caswell, W.E.: Asymptotic graph counting techniques in 
i/) 2Ar field theory, J. Math. Phys 19 (1978) 2579; Bender, CM., Cooper, F., 
Guralnik, G.S., Sharp, D.H., Roskies, R. and Silverstein, M.L.: Multilegged 
propagators in strong-coupling expansions, Phys. Rev. D 20 (1979) 1374. 



Breg.tex; 1/02/2008; 7:12; p. 9 



10 Solomon, Blasiak, Duchamp, Horzela, Penson 

11. Vasiliev, N.A.: Functional Methods in Quantum Field Theory and Statistical 
Physics, Gordon and Breach Publishers, Amsterdam, 1998. 

12. Blasiak, P., Penson, K.A. and Solomon, A. I.: Dobiriski-type relations and the 
log-normal distribution, J. Phys. A: Math. Gen. 36 (2003), L273. 



Breg.tex; 1/02/2008; 7:12; p. 10 



